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^ ■ 1 Statement of the problem and the main result 

Classical solvability of the Stefan problem for uniformly parabolic equations has 
O ■ been well studied - see for example papers [Ij - [9j and the references therein. 

At the same time, as has long been known, the heat transfer model based on 
uniformly parabolic equations, has some properties which can not be observed in 
, . , the reality, in particular, the infinite speed of propagation of disturbances. We 

I also know that more accurate model of the heat transfer is the model which is 

based on degenerate parabolic equations, such as equations of the form 

ut{x, t) = V{\ur-'Vu{x, t)) = fix, t), (1.1) 

where m > 1. As it is known, a short formulation of a classical Stefan problem 
for the equation (II. ip is the equation 

i(3{u))t = Vi\ur-'Vuix,t)) = 0, (1.2) 
where f3{u) is a discontinuous function of the form 
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u + k, u > 0, 



where > 0, is the latent heat of fusion (crystalhzation) and the equation fll.2p is 
considered in the sense of distributions. At that for a quasihnear equation (11. 2p 
the main unknown is, in fact, the interface {u = 0}, outside of which the solution 
of (11. 2p is smooth in view of the well-known local theory of uniformly parabolic 
equations. 

In its generalised formulation this problem was considered in a number of 
papers, from which we mention, for example, ^Oj - [16j . and we do not pretend 
to be complete in this matter, as the subject of our interest in this article is a 
smooth solution of the problem. 

As for the smooth solutions, in the case of one spatial variable such problem 
for degenerate parabolic equations was considered in [T7] - [21], where it was 
proved the existence of classical solutions. See also [22j . 

The aim of this paper is the proof of the classical solvability of the Stefan 
problem of the type (11.21) for a degenerate equation in a multidimensional setting, 
that is, the proof of the existence of smooth surface which is the interface {u = 0}, 
and the proof of the smoothness of the solution up to the interface. 

We now formulate a precise statement of the Stefan problem in a more ex- 
panded than (II. 2p form, as it is custom in the theory of free boundary problems. 
Let is a doubly connected domain in , whose boundary consists of two 
smooth connected surfaces and F" without self-intersections, dQ = F+ U F^. 
Suppose, further, that F is a given smooth surface without self-intersections ly- 
ing strictly between F"*" and F~ and separating the domain Q into two doubly 
connected subdomains and , so that dQ^ = F U F^. For a fixed T > 
denote = ^ x (0,T), n$ = x (0,T), Ft = F x [0,T], F^ = F± x [0,T]. 

Denote by St a smooth surface in the cylindrical domain Qt in the space 
(y, r) G X [0,T], such that at r = 0, it coincides with F, that is, S'^ fl {r = 
0} = F, St does not intersect surfaces F^ and divides the fi^ into two subdomains 
and Q^, and the lateral boundaries of these domains consist of 5*^ and F^ 
respectively. Surface St is unknown and has to be determined together with the 
functions u^(y,T) and u~{y,T), which are defined in respectively. The triple 
{St,u~^ ,u~) must satisfy the following conditions (we denote the independent 
variables by the {y, r) in view of the subsequent change of variables): 

-^-Vy{a^\u^r-'VyU^) = 0, {y,T)eQ$, (1.3) 
u+{y,T)=u-{y,T)=0, {y,T)eST, (1.4) 
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a+ Y,^os{f ,yi)\u+r-'u+ -a- J2^os{f, yi)\u-r-\^ = A;cos(i^,r), (y, r) e St, 

1=1 i=l 

(1.5) 

u^iy,T)=g^iy,T), iy,T)eT$, (1.6) 

M±(i/,0) = M±(|/). (1.7) 

Here m > 1, A; > 0, > are given constants, g^{y,r), g^{y,T), UQ{y), UQ{y) 
are given functions, at that 

±9^{y,r) >u>0, iy,r) E r±; ±4(y) > 0, y E <(i/) = 0, y G T, (1.8) 

where u is some positive constant: here and below we denoted by the same 
symbols u, /i, C all absolute constants, or constants that depend only on the 
given data of the problem. Note that the conditions (II. 4p . (ll.Sp are the three 
independent conditions at unknown boundary St arising from the equation (11. 2p . 

To formulate the smoothness conditions which we impose on the data of the 
problem, we introduce some weighted function spaces. First of all, we use the 
standard Ho"lder space W+^in) = C^+\n), 5 E {0,l),l E NU{0}, with the norm 
which are introduced in [23J, and also spaces H'''^^'^ (Qt) = C''^^'^ {^It) 
of functions of {y, r) with the norm In [23j the surface of the correspond- 

ing classes are also defined. We assume that the surfaces F, F^ belong to following 
classes 



F, F± E H^^^ (1.9) 

with some < 7 < 1. At the same time we suppose that the functions g"^ in 
(II. 6p are such that 

h^{y,T) ^ \gHy.r)r-'9Hy.r) E ff^+^'^(F±). (1.10) 

Suppose, further, that d^{y) is a given function from H'^'^"' {Q^) , which models 
the distance from a point y E to the surface F, that is 

< d^{y)/dist{y, F) < u-\ (1.11) 

Note, that such function can be taken, for example, as the solution of the following 
problem 

Ad+iy) = -1, yE n+, 
d+{y) = 0, y G F, d+{y) = I, y E F+. 
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Let, further, d (y) is an analogous function for the domain 
Denote here and below 

« = ^e(0,i), (1.12) 

m 

where m > 1 is the exponent from the equation fll.3p . 

2+7,— — ± rm 

We will use the spaces Cs ' ^ (^t) from the paper [H] (they are analogous 
to the corresponding spaces from |25]), where < 7 < a is some exponent, and 
we require 

0<7<a. (1.13) 

These spaces are defined in the following way. First we define the spaces 
C^s^^{R^t) in the domain 



R^^ = X [0,T], R^ = {x = {x',Xn) ■■ Xn >0,x' E R^'^} . (1.14) 
Define a distance between points x, x G R^ according to the following formula 

S{X,X)= /2 J/2 ' • (1-15) 

Define further a Holder constant of a function u{x, t) with respect to the variable 
X according to the distance we have introduced 



rr^ /x_ M -u{x,t) 
^sR^ M = sup -— . 1.16 



Denote by Cr^^iR^^) the space of functions u{x,t) with the finite norm 

|M|c.7,7/2(^iV ) = l^l^^^iv = \ufl^l +H]^j^M (u) + (M)J'^^iV , (1-17) 



(-) 

where (M)j^iv is the Holder constant with respect to t with the exponent 7/2 of 



the function u{x,t). 



Difine further the space ^ (R+t) ^s the Banach space of functions u{x, t) 

with the finite norm 



II I 1(2+7) I I I I 

'^'^ i=l 
N 

+ \Ut\c]''>/^^RN ) + \^NUx,Xj\c2'-f/\RN y (1-18) 



«J=1 
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Finally, the spaces are defined as the spaces 

of functions u{x,t) with the property, that in some neighborhood of Ft af- 
ter the corresponding change of variables functions u{x, t) belong to the space 

Cr^^R^r) or to the space Cs ^ (R+t) correspondingly, and out of some 
neighborhood of Tt the functions u{x, t) belong to the standard spaces C'^''^^'^{rij<) 

or C^^'*''^(f2^). In particular, for a function u G Cs'^'^' ^ (^t) ^^e following 
norm is finite 

\{d^)"u.^^^^\^^,-,/2^^±^ < oo, i,j = 1,N. (1.19) 

In the case of functions u{x) from the variable x only, x G the spaces 
C]{U ) and C2+T(n ) are defined in the completely analogous way. 

We will use also some standard anisotropic spaces of smooth functions, which 
are more general than the spaces if'''/^ = C'''/^. Namely, we will use the spaces 

where /i, I2 are noninteger positive numbers. Such spaces are defined in [26], for 
example. They consist from the functions, which have Holder continuous with 
respect to x withe the exponent li — [li] derivatives with respect to x up to the 
order [li] , and also these functions have Holder continuous with respect to t withe 
the exponent I2 — [I2] derivatives with respect to t up to the order [I2]. The norm 
of such space we denote by 



\U\^i,i„,-±, = \u 



ih,h) 



In fact, the functions from such spaces possess the property, that their deriva- 
tives with respect to x are smooth with respect to t and their derivatives with 
respect to t are smooth with respect to x. More precisely, if ki is a multiin- 
dex, that is, ki = A;i,2, ^i,Af) where ki^i are nonnegative integers, and 

k2 is a nonnegative integer, then the function 9^^(9f^M(x, t) belongs to the space 
C"^^''^^{Qj,) with the exponents rrii = /x/j, where 



/i = l- 



ll I2 



(see [26|), and also 



< C\u\^H,,^^.y (1-20) 

In addition, we will use the spaces with zero at the bottom of their notation, 
that is, the spaces (compare [23j, Ch.IV) 

clT'^iT)^), C[''\^i\ Cl'}{T)^l C',^'\Tt). (1.21) 
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Such notations means closed subspace of the corresponding space, which consists 
of functions that vanish at t = together with all their derivatives with respect 
to t up to the highest possible order in the corresponding space. 

Let here and throughout below the exponent /3 G (0, 1) is connected to the 
exponent 7 by the equality 

/5 = 7(l-f) (1-22) 



In particular, we will use the spaces C^^^~"'^"'"2 and C^^'^~"'^"'"2 (Ft), which 
consist of the functions t) with smoothness with respect to x up to the order 
2 + /3 — a and with smoothness with respect to t up to the order 1 + 7/2, that is 

N 
i=l 

It is easy to calculate, that because of the relation f ll.22p . we have 

7 2 + /3 - a , , 

1 + - = ^- , 1.24 

2 2-a ' ^ ^ 

so we also use the following notation for the mentioned above spaces 

C^+^--'^(nJ), C2+^-"'^(rT). (1.25) 
Further, we suppose that the initial conditions in (11. 7p are such that 

v^{y) ^ K{y)r-'u^{y) e Cl^^'m, (1.26) 

where 7' > 7. Besides, we suppose that 



dlt 



>^/>0, i/GT, (1.27) 



where 7^ is a normal vector to the surface F which is directed into 

We will show below, that the free (unknown) boundary St in fll.4p . (11. 5p can 
be parameterized in terms of its deviation from the given surface F^ = F x [0, T]. 
We follow to [2j to give the strict formulation. Let uj = {ui, ...,Un-i) is a local 
curvilinear coordinates in a domain B on F. In some small enough neighbourhood 
TV in of the surface F we introduce the coordinates {u, A) in the way that for 
any x G A/" we have the following unique representation 

X = x'{x) + lt{x'{x))X = x{u) + lt{u)X, (1.28) 

where x'{x) = x{u) is the point in the domain on the surface F with the 
coordinates u, lt{ijj) - normal to F at the point xiuj) with the direction into Vf^ , 
A G i? means, in fact, deviation of the point x from F, at that ±A > for x G VL"^. 
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We assume that the mentioned above neighbourhood of the surface F is the 
set 



J\f={xen: \X{x)\ < 70} 



where 70 is small enough and will be chosen below. 

Let p{x',t) = p{u},t) is a small and regular function, which is defined on the 
surface Tj'. Let us note, that here and in what follows we use the notation p{uj, t) 
with the argument u instead of p{x', t) for all functions on the surface F if it does 
not cause ambiguity. We do that just for simplification of the notation, bearing 
in mind that in each local domain B on F we can introduce local coordinates u. 
At the same time the coordinate A in flL28p does not depend on the choice of 
local coordinates u. 

We parameterize the unknown surface St we the help of the unknown function 
p{u, t) as follows 

St = Tp^T = {ix,t) G fir : x = x' + p{x' ,t)lt{x') = x{u)+ p{u,t)lt{u)}, (L29) 

where x' = x{uj) G F. Note, that this definition of the surface St = ^p,T does 
not depend on a choice of local coordinates u in a particular local domain on 
F. Thus, the unknown function p{u},t) means, in fact, deviation of the surface 
Tp.T = St from the given surface F^. 

Along with Q^, in (II. 3p we use the notation = and ^~t = Qt f*-"^ 
the subdomains that Tp^T = St separates the domain Qt- Let, further, p{x,t) is 
an extension of the function p{u, t) from the surface F^ to the whole domain Qt 
to a function with the support in the neighborhood At = A/'x [0, T] of the surface 
Tt, p{x,t) = Ep{u,t), E is some fixed extension operator (the way of such an 
extension will be listed below), at that we will denote p"^ = Ep\-^± = E^p. 

Define a mapping ep{x,t) from the domain Qt on itself with the help of the 
formula Cp : {x,t) — {y,T), where, according to the notations of f ll.28p . 



Here x'{x) G F, uj{x), X{x) are (w. A)- coordinates of a point x in the neighbour- 
hood JV. Note here, that the definition of the mapping Cp does not depend on a 
choice of local coordinates u on the surface F. 




(L30) 




x'{u{x)) + (a;(x))(A(x) + p{x, t)), x e A/", 
X, X eQ\Af, 

T = t. 



(1-31) 
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We choose 70 small enough so that under the condition 

Iplr^' < 270 (1.32) 

the mapping Cp is a diffeomorphism of fij- on themselves, and also the mapping 
Cp is a diffeomorphism of the domains Vt^ on the domains i^py. Note, that the 
surface Fp 7^ is exactly the image of the surface Tt under this mapping, and the 
mapping ep{x,t) is the identical mapping out of the neighbourhood A/j- of Ft- 
About the exponents of the Holder spaces we use we suppose that 

< 7 < 7' < 1, 7 < min{Q;, 1 - a}. (1.33) 

Note, that under our chice of 7 the restriction 7 < is also fulfilled. The 

last restriction was introduced in [21] at the studying of the homogeneous Cauchy- 
Dirichlet problem for degenerate equations. We need the restriction 7 < 1 — a to 
have the inequality 1 + — q;>1 — q;>7, which implies that the first derivatives 

2 13 a ^+^~" 2+7, — -I- 

with respect to x of the functions from the classes C ' , Cs ' ^ (^t) 
are more smooth than C'^''^/'^^ C]'"'^'^. At last, we need the restriction 7 < a 
to expressions of the form {d^{x))'^ri{x,t) with a smooth function rj{x,t) would 
belong to the corresponding space C]'"'^'^. 
Let us formulate now the main result. 

Theorem 1.1 Let the conditions f ll.8p -( TTTTU]) . (I1.33P on the data of the problem 
(11. 3p - (11. 7p and the conditions (I1.26P . (I1.27P are fulfilled. Let also the natural 
adjustment conditions up to the first order at r = 0, y G F, for the problem 
(II. 3p - (ll.7p are fulfilled. Then there exists such T > 0, that on the time interval 
[0, T] the problem (II. 3p - (Il.Tp has the unique smooth solution, at that the the 
unknown boundary can be represented as in (I1.29P with the function p{u, t) with 
the properties 

p{u,t) G C2+^-"''^^(Ft), pt{u;,t) G C1+^-"''^^(Ft), (1.34) 

p±(x,t)=E±p(a;,t)GCf^(nJ), 

where p{x,t) = Ep{uj,t) is the extension of the function p{uj,t) to the domain 
Qt- The functions u'^{y,T) in (II. 3p are such that 

v^{x,t) ^ {\u^{y,T)r-\^{y,T))oep{x,t) G C'+'^iHt). (1-35) 

Thus, in particular, all of the relations of the problem (II. 3p - (ll.7p are satisfied 
in the classical sense. 

Subsequent sections of the paper devoted to the proof of the theorem 11.11 in 
accordance with such a plan. 
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In the section [2] on the basis of equivalent norms in the spaces 

the natural space of the traces on of the functions from class Cs ' ^ (^t) 
is studied. This allows us to extend the results of [21] about the solvability of 
the homogeneous initial boundary value problem for a degenerate equation to the 
case of the inhomogeneous problem. 

These results are then used in the section [3] to study a model Stefan problem 
for degenerate equations, which is one of the central points of this paper. In this 
case, for the Schauder estimates of the model problems the idea of the paper |28] 
on the application of the maximum principle to obtain the Schauder estimates is 
used. 

In the sections H] and |5] the initial problem with the unknown boundary is 
reduced to nonlinear problem in the fixed domain and linearized on functions 
that expand the initial data to the domain t > 0. 

The section |6] is devoted to the study of the resulting linear problem for de- 
generate equations on the basis of the results of the section [3] about the properties 
of the corresponding model problem. In this case, to prove the solvability of the 
linear problem, we apply the idea of [1] on parabolic regularization of the Stefan 
boundary condition. Note that the corresponding model problem in Section [3] is 
considered in the presence of the regularization. 

Finally, the section [7] completes the proof of the theorem 11.11 by the method 
of i3j. 

2 Auxiliary results on the spaces Cj'^'^ (^Qji), Cg {Vtj^ J ^ 

2+/3-a 



Note first, that for the spaces with zero in fll.2ip the following relations are valid. 
Let 7' > 7, l'^ > li, I2 > h- Let also functions u and v belong to one of the 
mentioned spaces with the exponents of smoothness 7', 2 + 7', l[, Zg- Then 

\uf"'l^ < CT^\4^2^, (2.1) 
|^|(2+7) < cT^\^f'l2'\ (2.2) 

S J ^^rp S , ^^rp 

\uvf''l^<CT^\uft^\vf"'l^, (2.3) 

^^^{h,h) ^ CTf'\ufY'\ (2.4) 
where fi is some positive constants which depend on 7, 7', li, l[. 



These inequalities are well known for the spaces Cg^' ^ (see [23], [27]), and for 
the spaces Cq']^"^, C^'^.'^' ^ these inequalities are completely analogous. 
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2.1 An equivalent norm for the spaces C]{^^^), Cg~^^{Q^). 

Along with the seminorm -f^J^jv from f ll.lGp we consider in R^j^ the following 
weighted seminorm 



Hla^ (/)= -P '^'Y'"-|f'' " + 4* '^'"-,""4f'" . (2.5) 



where x^- = max{x7v, a^Af}, and here and throughout we, without loss of general- 
ity, assume that Xat < x^, so x^ = x^. 

Lemma 2.1 The seminorms j^n (/) and Hj (/) are equivalent. 
Proof. Let the seminorm t^n (/) is finite. We show, that then 



+T 



KR^Jf)<CHl^^Jf)- (2.6) 
Let ^ (0, 1) is small and fixed. Let x = {x',xj\f) and let first 

\x' — x'\> SqXn- (2.7) 

Then the more 

|x — x| > \x' — x'\ > SqXm- (2.8) 

Under this condition 

s(yX^ all '~~ 

X' I />'■ ' I I <-y»/ rv*^ \ OL I '2i 

j\j ~[~ ]\[ 1"^ I 

3/ ^ I ^ ^ I I -I a. 



< C ' — < C\x-x\'-^ 



Therefore, as /3 = 7(1 " 



2 ^' 



'^(^'^^ < < CH\^ if). (2.9) 

|x — x|p s{x,xp *'-"+T 

Besides, because of fl2.8p . and then of (12. 9p . 



^,^\f{x,t)-f{x,t)\ ^ xl^ \f{x,t)-f{x,t)\ 



^^^^^ '-^-^ < '\ . 'V ' \ n /^V < CH\^ if). (2.10) 



Let now 



10 



\x' — x'\ < EqXn- (2-11) 

Under this condition, as it easy to see, 

s(x, x) ~ CxjY^ |x — x|. (2-12) 

Consequently, 

^n \f{^ t)-mt)\ ^ ^ \fix,t)-mt)\ < (2.13) 

^ \x-x\^ ~ s{x,x)^ ~ ^.R^T^''^ ^ ^ 



To estimate, further, the unweighted Holder constant in the definition of 
H"' ff (/), we consider the two cases. 

Oi ^ JL j rj-i 

If 

\xn -xn\> £oXn, 

then 

|x — x| > |Xiv — > SqXjsi 

and therefore, as it was above. 



X — x\ , _ 



six.x) < — — < C\x — X 



2 



SO that, as above. 



< C '^^ ' " < CH:^^ if). (2.14) 



If now, under the condition fl2.1ip . we have 



\xn — xn\ < £oXn, (2-15) 



then in this case 



|x — a;| < \x' — x'\ + \xn — x^l < "^eoXN. (2-16) 
Therefore, in the force of ( I2.12p . 



s{x, x) < Cxj^"^'^\x — x\ < 



< Cx-''^\2eoXNT^^\x - = C\x - xj^-^/l (2.17) 

Consequently, in this case 
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l^'---'>-^f-"l < cl^'^-'>-{'^-'>l < CH^„. if). (2.18) 

The estimate (^B) follows now from 021]), fETTU]) . (I^TT^ . (l^Ti]) and d^TTg]) . 

Further, let now the seminorm H^j^n if) is finite. Let us prove the following 
estimate 



Let first 



Then 



and consequently 



KR^Jf)<CHln!tM (2-19) 



\x' - x'\ < EqXn, xn>0. (2.20) 



s(x, x) > V- 



\x — x\ 



a/2 



\f{x,t) - f{x,t)\ ^ ^^^a/2 \f{x,t) - f{x,t)\ ^ 



In the particular case a; at = we have x^ = and therefore 



o(rr ™N _ I ' _ 77T'|l-a/2 _ I _ 7-t|1-«/2 



and so again 



s(x, |x — X 

Let now we have 



-+T 



\x'- 


x'\ > 


\x ■ 


- x\ 


\x' — 





(2.21) 



^1, <CHln^Af)- (2.22) 



(2.23) 



Then 

six,x) > v ^J j;,, > z/|x-a;|^-"/^ (2.24) 
and consequently, 

'^(-^■''-{'^■"l < c^f<^A-f^A\ < c„^ ^„ (2.25) 

Thus, flCTj) follows from ^TH) . ([222]), (ESI. And so the equivalence of the 
seminorms H^^ if) and H'^ (/) is proved. ■ 
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In this way, the norm in the space C]{R^rp) may be given in the form 

= l^lSi + Kn^r^"^^ + (^)S^^- (2.26) 

Bearing in mind the local straightening of the boundary F, for the case of 
arbitrary domains fi^ , the norm in the space may be explicitly written 

as 



(2.27) 

where the functions d^{x) were introduced in the previous section in fll.lip and 

they model the distance to the boundary F, d^{x,x) = ma:x{d^{x),d^{x)}. 

Quite similar in terms of f|2.27p and ( I1.19P we may explicitly define the norm 
ffLTSjl in the space C2+T(nr). 



2+7 ^ ± 

2.2 The traces of the functions from Cs ^ (^t) ^t- 



In view of the smoothness of the surface F, we can use local straightening of 

2+7, — ± 

the surface at consideration locally defined classes Cs ' ^ (fiy). Therefore it is 
sufficient to consider the case of the half-space, that is to consider the finite in 

2+7 

-R+T = -R+ ^ [O5 ^] function u{x,t) from the space Cs ' ^ {R+t) consider 
its trace at a; at = 0. 



Lemma 2.2 Let the function u{x, t) IS finite and u{x,t) G C^+^'^{R^t) , < 



^ <:a, ^ = 7(l-a/2). Then the function v{x\t) = u{x',0,t) G 



2+/3~-a, ^+g'° / pjV-l^ 



C2+/3-,^(i?^-i)^ at that 



T 



Hx\t)\ = \uix',0,t)\ <C\<V ■ (2-28) 

C ^ ' (fly ) C ^ ' (fly ) s,n_,_y 

Besides, 

|V(.,.,)«(x',0,t)|^,,,_i+^^^^_^^ < (2.29) 

Proof. It follows directly from the definition of the space 
ffLTSj) and from the lemma EIH that Ut{x',0,t) e ^^•'^/^(i?^"^), and therefore 
vtix',0) = utix',0,t) e C'^'^/'^iR^-^), and in addition 

l^tlc^^V2(ij,^-i) < C \u\f^^^^^ . (2.30) 
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Therefore, in the force of (ll.20p (see [20]), it is sufficient to prove uniformly in t 
the following estimate 



'■+T 



(2.31) 



and for this it is sufficient to prove, that uniformly in t and in Xat for alH = 1, 
for the function w = u^. we have 



, ai+/3-a) < r.L,|(2+7) 



+T 



(2.32) 



So, let w = Ux^, i = 1,N. To prove f l2.32p it is sufficient, as it follows from 
to show that for arbitrary h > the follows inequality holds 



(2.33) 



Here A\ ^,w{xn) = '^Ix''^ = ^\x''^i^' i^N,t) is the second difference from the 
function w with respect to the variable x^-, j = 1, — 1, with the step /i, that is 

^iw{xn) = u{x' + ~j h,XN, t) — 2u{x', XN,t) + u{x' — ~] xn, t). (2.34) 
Consider the two cases. Let first 



h < xn- (2.35) 
Then according to the mean value theorem with some 61,62 G (0, 1) we have 

\AIx,w{xn)\ ^ \x%AI^,w{xn)\ _ 



X 



Wx, {x' + 7" 6^1/1, XN,t) - Wx, {x' - ~j'62h, Xn, t) 



N' 



h0 



x^Ux^x, (x' + ~^6ih, Xn, t) - x%Ux^xi {x' + ~t02h, xn, t) 



< 



< C {x%Ux^x,)^^}^^N^ < C • 
Let now h > Xn- Write the difference Aj^ ^,w{xn) in the form 

Al ^,w{xn) = -Al ^, {w{x' , Xn + h,t) - w{x', xn + h, t)) + 



(2.36) 



+ Al x'W{x', Xn + h,t) = Ai + A2. 



(2.37) 
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In view of the fact that for the expression A2 the condition fl2.35p holds, that is 
h < xj\f + h, completely analogous to (12.36 p . 



<C\u 



(2+7) 



(2.38) 



To estimate the expression " we use the formula 



w{x' ,Xn + h,t) — w{x' ,Xiy,t) = h J Wxr^{x' ,Xiy + 9h,t)d9. 





Consequently, in view of w 



< 







Al ^, {xn + OhYu^^^^ (x', xn + 0h, t) 



dd < 



< C {x%u,,.Jf,]^r^^ I {x^/h + e)-''d9<C . 



(2.39) 







Thus, from fl2.38p and f l2.39p we obtain fl2.33p . and so we have also fl2.32p and 
((231]). Together with ^LM) this completes the proof of (E28]). 

We now show the inequality fl2.29p . Note that for tangential derivatives Ux^, 



i = l^N — 1 this inequality follows from the above estimate fl2.28p and from 
f ll.20p . However, for u^j^ we need a separate proof. We show fr2.29p for m^.^.. 

According to [29j, it is enough to show that for h > 



\Al,Ux^{x,t)\<C\4^^;i>h 



1+13- 



+T 



(2.40) 



where 



Al ^v{x, t) = Al^v = v{x, t + 2h) - 2v{x, t + h) + v{x, t) 



is the second difference of a function v with respect to the variable t with the 
step h. 
Let first 



/i2-" < Xn- 

Then we use the following interpolation inequality (see, for example, 
Ch.l) 



(2.41) 



\^^k\nT(xM) - ^ U^lnT(xiv) 



1/2 



I (2) 



1/2 
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Here we denote IIt{xn) = {{y,t) : Xn/2 < i/n < '^Xn/'^^^O < t < T}, and |f 1^^,(3,^) 
means C^- norm with respect to x- variables over the specified domain. We obtain 
for Al^U:rk{x,t), that 

|A^,«.,(x,t)| < C {lAlufZ^T' (2.42) 

In view of the properties of the space Cs'^'^' ^ {R+t) (c^- EH], [21], (H.20P ). we 
have 

|AL"l£,«, < CWClh^. (2.43) 

< C\u\<^_;lx7rhi. (2.44) 
From fl2.43p and fl2.44p considering f l2.4ip . we obtain 



'■+T 



(2.45) 



that is the inequahty (12.401) . 
Let now 



/i2-« > Xn. 



Write ^u^^^ (x, t) as 



ALwxfc (a;, t) = - AL {x\ xn + h 2— , t) - u^^ (x, t) 



+ 



(2.46) 



+ Al^u^^ {x\ Xn + h^-'-,t) = Ai + A2, 
and for A2, by the above case (I2.4ip . the estimate 



is vahd. 

To estimate the expression Ai, write it as 



Ai = -h^ j Al^u^^^^ (x', Xn + Oh^ , t)d9. 



Thus, we have for Ai, that 



(2.47) 



\Ai\ < C\u\^^p^h^hi I (xn + 9h^ ) d9 < 
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S ^ 1\ j rj-1 J S ^ jTh j 



5 



that is again the inequahty (I2.40p . 
The lemma is proved. 
■ 

Thus, due to the possibihty of the local straightening of the boundary, the 
following is true. 

2+7,2+2 —± 

Lemma 2.3 Let functions u {x,t) belong to the spaces Cs ' ^ (fi^). Then the 
functions v'^^Xyt) = u^{x,t)\x£r belong to the space C^^^^°' ^-q (Ft), and 



\v^\ 2+^- <C\ufl2^. (2.48) 



Besides, 



|V«^|r, I ^ < C\u^\%2^. (2.49) 

C ^-O! (r^) SAlrp 

2.3 An extension from the surface Ft of the functions from 
the space C^+^-'''-^ {Ft). 



In this section we prove the converse of Lemma [273| that is, we show that any func- 
tion of the class 0^+^-""'^^ [Ft) = C^+^-'^^^+^/^Tt) can be extended to all do- 

-I- 2+7,2±21 ± 

mains Qrp can be extended to all region up to functions of the class Cs ' ^ (^t); 
and the extension operator is bounded (here, as above /3 = 7(1 — a/2)). Such 
an extension operator is constructed in the standard way by applying a suffi- 
ciently small partition of the unity in the neighborhood of F and by the local 
straightening of the boundary F - see |23| . In this case, it is enough to require 
the iJ^+^-smoothness of the boundary F. Therefore, the existence of the said 
extension operator follows in the standard way from the following lemma. 

Lemma 2.4 Let in R'^j, at = a finite function f{x',t) from the class 
Q2+(3-a,i+y/2^j^N-i-^ zs givcn. Thcn f can be extended in the domain Xn > 

up to the function u{x,t) from the class Cs ' ^ {R^rp), and 

< C . (2.50) 

Proof. Let u{x,t) is the solution of the following Dirichlet problem with the 
parameter t E [0,T]: 

Am = 0, xeR^ {xn > 0), (2.51) 
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uU^=o = fix',t), (2.52) 

u-> 0, \x\ oo. (2.53) 

As it is well known, the solution of (I2.5ip - f l2.53p is given by the potential of the 
double layer, defined by the Newton potential. 

Note, first, that for the problem (12.511) - (I2.53P we have the following maximum 
principle 

< l/lgL.. (2.54) 

Indeed, by (I2.53p . we can choose K > so large that \u\ < |/|^|r_i/2 for |a;| > K, 
and, by the properties of the double layer potential and the finiteness of f, can 
be chosen independent of t. Now consider in the domain = -R+ {1^1 < K} 
the functions v"^ = ±u + |/|^lr_i- Within this domain we have 

Av^ = 0, X e Bk, (2.55) 

and on the boundary dBx = {xn = 0, < K} U {xn > 0, \x'\ = K} the in- 
equality 

v^>0, xE OBk (2.56) 

holds. It follows from (I2.55p . (I2.56P and from the maximum principle, that f ^ > 
in Bk for all t, and, thus, 

Hgi<l/l2L. tG[0,T]. 

Due to the choice of K, we have the inequality (12.541) on whole domain R^rp. 

Consider first the properties of the function u{x,t) with respect to t. Let 
v{x,t) is the solution of fl2.5ip - fl2.53p with the boundary condition 

v\,,=o = Mx\t) (2.57) 
instead of f l2.52p . Consider also for /i > the following function 

u(x, t + h) — u(x, t) 
Uh{x,t) = , 

which satisfies the problem fl2.5ip - f l2.53p with the boundary condition 

, . / ,^ _ fix',t + h)-fix\t) 

Uh\xN=o = Jh[x , t) = . (2.58) 

Let further w{x,t) = Uh{x,t) — v{x,t), and the function w{x,t) also satisfies the 
problem fl2.5ip - f l2.53p . but with the boundary condition 
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wUr,=o = fhix',t) - ft{x',t) = iph{x',t). (2.59) 

Due to the properties of the function f{x\ t), we have with some 6{x', t, h) G (0, 1) 
according to the mean value theorem 

\(phix',t)\ = \fhix',t) - ftix',t) \ = 



\ft{x',t + eh) - ft{x',t)\ < {ft{x',t))'^;^l^l, W''' ^0, h-^Q. (2.60) 



Consequently, on the base of f l2.54p . 



(2.61) 



This means, that the function u{x,t) has the derivative with respect to the vari- 
able t for X G -R+j., and Ut{x,t) = v{x,t), that is Ut{x,t) satisfies the problem 
fITM]) - 0233]) with the boundary condition ([22ZD- 
Further, considering the function 



Vh{x,t) 



Ut{x, t + h) — Ut{x, t) 



we see, that it satisfies the same problem with the boundary condition 

ft{x',t + h)- ft{x',t) 



Vh{x,t%j^=Q 

Thus, on the base of (I2.54p again. 



fth{x',t). 



(2.62) 



Ut{x, t + h) — Ut{x, t) 



(0) 



< C 



ft{x',t + h)- ft{x',t) 



.N-1 
T 



which, by the arbitrariness of /i, means that 



(2.63) 



(2.64) 



{<^MSl<C{ft{x\t))^[l^l,. 

Consider now the properties of the function u{x,t) with respect to the vari- 
ables X. 

First, it follows from the results of [21], that for each t G [0,T], due to 
/ G C2+/^"°(i?^-i), we have u G C^+''-"'{R'^t), and 



(2.65) 



where the symbol x at the bottom of the space notation means that we consider 
the smoothness only with respect to x. 
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Show that the following estimates 

{^ND\)i%^^ < (2.66) 

(xl^/' {^%D\)f^^^^^^ < (2.67) 

are valid, that is 

H2{x%D^u) < C|/|^2+,-.(^^-i). (2.68) 

We will use the fact that, as it follows from [33], Ch.5.4, the condition / G 
Cl{R^-'^) in fl232|) . / G (0,2), is equivalent to the condition 

\D'M < CkX],'^'\f\cUR^-i), k > 2, (2.69) 

where here and below D^u = D^u means a derivative of the fc-th order with 
respect to x of the function u{x,t). 

Since it is important to prove (12.661) for xat < 1 only (for xat > 1 such the 
estimate follows from the local estimates and is well- known), we consider only 
the case x^ < 1. 

We also use the well-known interpolation inequality 

("M>Ssc(MS')'"(l4')''. ("«) 

which is valid for the functions v{x) G C^{Q), Qis a (possibly unbounded) domain 
with the sufficiently smooth boundary (see, for example, [30], Ch.l ). 

In addition, at the proof of (12. 66 p . without loss of generality, we prove smooth- 
ness of the function x'^D'^u with respect to x' and with respect to Xn separately 
and we obtain the estimate fl2.66p separately for these two cases. 

So, let first Xn is fixed. Then, by flZTO]) and flT^ . 

{x%D\)l',^ < C {\x%D\\('^Y'' {\x%D\\'^'^ + \V..x%D\rY < 



(2.71) 

that is estimate (I2.66P with respect to x'. 

Analogously, using fl2.70p and f l2.69p . we prove f l2.67p with respect to x': 

xl^/' {x%D\)^) < Cxl^'^ {\x%D\rf~' [\x%D\^y < 
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^ r\ fl /^^" ^-2+(2+/3-a)y / „ 3+(2+/3_a)y _ 

= C\f\^2+P-<.^j^N-l^x]^^'^^^ ^ = C\f\^2+P-c^j^N^ly (2.72) 

We prove now the relations (I2.66p . (I2.67P with respect to the variable x^. For 
this we fix some eo G (0, 1/16) and consider the two cases, assuming without loss 
of generality that xn < Xat. 

Let first 

\xn — xn\ = — xn) > i^qXm- (2.73) 

Then 

\x%D'^u{x, t) - x%D^u{x, t) I 



\Xn — Xj\f 

In this case, as above 



/3 



< C [\x%-''D\{x,t)\ + \x%-^D\{x,t)\^ . (2.74) 



x% ^D'^u{x,t)\ < C|/|^2+/3-c(^iv-i)X^ ^a;^^^^^"^'^ = C\f\(j2+i3-c^j^N-iy (2.75) 



and similarly for ^ D'^u{x,t)\. 
In the same way 



7a/2 



x%D^u{x, t) - x%D^u{x, t) I 
\xn - XnV' 



and then proceeding as in f l2.75p . 
Let now 



< C (\x''-^D\{x,t)\ + \l!f'-^D\{x,t)\) 

(2.76) 



< {xn — xn) < £qXn, (2.77) 

and let also 

n(xjv) = [y e R1 : xn - 2eoXN <yN <xn + S^oa^jv} , 

nT(x,v) = n(x^) X [0,T]. (2.78) 

Then, taking into account that on IIt{,xisi) we have ~ a;Af, as in the previous 

case 
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< C [\y''^-'D\tZ^^^^ + x% {DMy,t))l%^^^^^) ^A, + A,. (2.79) 

Here Ai is estimated in the same way as in (I2.75p . and A2 - as well as in f l2.7ip . 
which gives 

\x%DMx,t) -x%DM'3^,t)\ 



Xjv — X]\i I 



The estimate 



■ya/2\x%D'^u{x,t) X'^D'^u{x , t)\ ^ „: /o CI ^ 

is quite similar. 

This completes the proof of f l2.68p . 

Similarly, we obtain the properties with respect to the variables x of the 
derivative Ut, that is, 

H2{ut)<C\f\^2^,^.^^M-.^ (2.82) 

because Ut satisfies the problem fl2.5ip - f l2.53p with the boundary condition (I2.57p . 
Indeed, since Ut\xj^=o = ft, so 

< C {ft)i%.. < C\f\^,,,_^._^^^^^_^. (2.83) 

Further, for x,x E R^, Xn > Xn consider the difference 

Aix,x)u, = ^.^^/ M^f)-Mx,t)\ _ ^2.84) 



If \x — x\ > EqXn, then 



< C\f\ 2+,8-. „ , . (2.85) 

If now \x — x\ < SqXni then Xat ~ Xat, and then, using fl2.69p . we obtain that 
A(x,x)«, < ^^W2 l^t(x,t)-y,t)| _ ^ 

X — X 
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(2.86) 

Now ^Mj follows from flXSHj) and (IZMj) . 

Let us show now the smoothness of the function x'^D'^u{x,t) with respect to 
the variable t, that is show that 



{x%D\)^'^^^^ < C|/|^,^^__2+^^^^_^^ = C\f\c2+p.^,i+-,/2(^RM-iy (2.87) 
For this we fix some h > and consider the function 



u(x,t + h) — u(x,t) 
t) = ^ J/2 ^ (2-88) 

which satisfies the problem 02.5ip - (I2.53P with the following boundary condition 

M / N f(x,t + h)—f(x,t) , , 

Vh{x,t)\,,=o = M^,t) ^ ^ • (2-89) 

It follows from the results of (SB], (ll.20p . that uniformly with respect to the 
variable t the function iph{x,t) G C^~"(i?^^^) with respect to the variables x, 
and 

max \(ph{-,t)f^^"l < C\f\c2+p-c.A+^/2(^RN-iy (2.90) 

Note now, that 

x%D'^u{x,t + h) - x%D'^u{x,t) ^ 2 / ,N 
= x^D Vh{x,t). 

Consequently, it follows from (I2.69p that 



\x%D'^Vh{x,t)\ < Cx% max \iph{-,t)fL'^l x,^^^^"^ = C max \(ph{-,t)fL'^l . 

t6[0,T] " te[o,T] 

(2.91) 

So, (IZSTj) follows from (ICTD and I^M^, in view of the definition of v{x,t). 

Multiplying now the function u{x,t) by a smooth finite function r]{x), which 
is equal to one in a neighborhood of support of f{x', t), we get a finite extension 
of f{x',t) of desired class, and the estimate f l2.5Up . 

The lemma 12.41 is proved. ■ 

From this lemma in the standard way (see |i23j). as it was described in the 
beginning of this section, we get the following assertion. 
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Lemma 2.5 There exist bounded extension operators and E , such that 



e c2+^-'^(rT) ^Ep^ E^p = p± e ^^^''^^(nj), (2.92) 



p 



and we can assume that the supports of the extended functions p^ are included 
in the sufficiently small neighbourhood of the surface Ft- We will denote the 
extended functions p = Ep = E^p by the same symbol p to not to overload the 
notation, that is, 

P\n^^P^^E^P^Ep\n^- (2-94) 

Besides, as it follows from the results of [21] and from the lemma 12.41 (as 
the lemma 12.41 permits to reduce the situation to the homogeneous boundary 
conditions), the following assertion is valid. 

Lemma 2.6 Let functions f and g are finite, and 

f{x',t) G ' g{x,t) G CI,{R^t). 

Then the problem 

— -x%Au = gix,t), ix,t)eR^T. (2.95) 

u{x', 0, t) = f{x', t), XN = 0,te [0, T], (2.96) 

m(x,0)=0, xeR+ (2.97) 
has the unique solution u{x,t), which satisfies the estimate 

|.|S?<c(|/|^„,_.,^^^„_.^ + Wi:i»^). (2.98) 

In the same way, with the help of results of [24] and from the lemma 12.41 we 
get the following theorem. 

Consider the Cauchy-Dirichlet problem for the degenerate equations of the 
form 

^ - {d^{x)yB^{x, t)Au^ = g^{x, t), (x, t) G (2.99) 
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u^\r^ = f^{x,t), (2.100) 

u\± = h^{x,t), (2.101) 

M^(x,0) = 0, (2.102) 
where the functions d^{x) are introduced in (II. lip . 

B^{x, t) e C^'^/^Qt), i'<B^ <u-\ 

G Cl-J\Tt^\ e cl^'~''^'^{Tr), G cl"-''"^ {T$). (2.103) 

Theorem 2.7 The problem (I2.99P - (I2.102P has the unique solution from the 
space Cq^ ' ^ (^r) '^'^^ following estimate is valid 

3 The model problem for the two phase Stefan 
problem for the degenerate equations. 

Let a > is a fixed number. Denote = {{x,t) : x G R±,t > —a}, Q^^^ = 
{{x',t) : x' G R^-\t > -a}. Denote further = Q± n {t < T}, R^'^''' = 

qN-i n {t < T}. It is convenient to consider the domains with the t > —a, as 
it will allow us to consider the points with t = as interior points of general 
position, which will facilitate the further notation. We agree, which is similar 
to ( ll.2ip . that zero at the bottom of the designation of the spaces of functions 
defined in these domains means the subspace of the corresponding space whose 
elements vanish at t = —a together with all its derivatives with respect to t, 
which are permitted by the space. 

Let f{x', t) is a finite with respect to x function, which is defined in Q^^^ and 
is such that 

/(x',-a)=0, /gCo'^''""''^(Q^"'), (3.1) 

which allows us to consider / as the functions, which is defined for t G (— oo, oo), 
extending it by identical zero in the domain t < —a with the preservation of the 
class. 

Let further 



(2.104) 
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fHx, t) e Cl^'^Qll ft{x\ t) e cI^'-''^^'\q''-') (3.2) 

- are given finite functions wliicli are also extended by identical zero in the domain 
t < —a with the preservation of the classes. 

Consider the following model problem for the triple of the unknown functions 
u^{x,t) and p{x',t), which are defined in Q_^_ and Q correspondingly: 

^-{±x^yAu^ = ft, {x,t)eQl (3.3) 

M± + A^p = f^, xn = 0, {x', t) G Q^-\ (3.4) 



(/l I (11 1 

p,-eA,,p + b+- b-^ = f{x\t), x^ = 0,(x',t)Gg^-\ (3.5) 



dx 



N 



dx 



N 



u {x, —a) = 0, p{x', —a) = 0, 



(3.6) 



u e C,_o (<5±), P e Co " [Q ),pteCQ [Q ), (3.7) 

where A^, b^, e are given positive constants and e G (0, 1). 

Note that the term with e in fl3.5p does not apply directly to the Stefan 
problem and serves as a regularization of the problem, that will be needed in 
the proof of the solvability of the corresponding linearized Stefan problem in an 
arbitrary domain. To the author's knowledge, this regularization of the boundary 
condition in the Stefan problem was first used in the paper pLj. 

Below we prove the following a priori estimate of the solution of the problem 

(E31) - (1X7|). 

Theorem 3.1 Let u^{x,t) G c'o^'^(Q±), p G C^^^'"'^^^^ (Q^~') are a 
finite solution of the problem ( 13. 3p - (13. 7p . Then for arbitrary T > the following 
estimate is valid: 



U{T) ^ + \u 



I (2+7) 



' I \^ I 



+ |pLiV-i,a 
I (7) 



+ \Pt 



R: 



N-l,a 



< 



,R 



+T 
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where the constant Ct in f l3.8p (ioes not depend on e E (0, 1). 



Subsequent content of this section is the proof of the theorem 13.11 
Note first that by lemmas [2 . 6 1 and we can without loss of generality assume 
that 

/±^0, ft^O, (3.9) 

since the general case can be reduced to the specified one by the change of the 
unknown functions = + w^, where are the new unknowns, and 
satisfy f l3.3p with the boundary conditions 

w\,=o = fHx',t). (3.10) 

Thus, further we assume that only the function f{x',t) is nonzero in the 
righthand sides of (13. 3p - (13. 5p . 

In addition, because the right side of the relations ( 13. 3 p - (13. 5p belong to the 
classes with zero at the bottom and because of conditions (13. 6p . ( 13. 7p we can 
consider that the relation (13.30 - (13. 5p are valid and for t < —a, assuming that 
all the functions are extended by zero to this domain. 

An important point of proving (13.81) is to prove the following a priori estimate. 



Lemma 3.2 Under the conditions of the theorem \ 3.1\ and under the condition 
(I3.9p the following estimate is valid 

= Ct {Af{T) + M{T)) < Ct{T + a)'^W(T) + CtM{T), (3.11) 



where 



To obtain the last inequality in (13. lip we use the estimates (12. ip - (12. 4p . 
Proof. 

Denote for brevity, I = 1 + f3 — a and fix a point (xo,to) in the set R^'^'"". 
In order to maintain the succession of the notations with the paper [28], whose 
method we're going to apply, without loss of generality, we will assume that 
(xo,to) = (0,0) - this choice is not important, as can be seen from the following 
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proof. Suppose, further, O = {x' = Q^x^ = 0,/!: = 0) is the corresponding point 
in R^^. We show that for every h G (0, 1) and for any j = 1, — 1 we have 
the following inequality 

\p.A'^,h, 0) - p.,(-^,/i, 0)1 < Ct (A/'(T) + M{T)) h\ (3.12) 

where ~^ j is the unit vector of the Oxj- axis. Since the point O and the step 
h G (0, 1) in the relation (I3.12p are arbitrary, the estimate (13. lip of the lemma 
follows from the estimate (I3.12p . 

So, let 2/1, 2/2 G [0, 1], ?/ = (yi, 2/2) and let also z, j G {1, 2, — 1} are fixed. 
Consider the differences 

v^{x,t,yi,y2) = A,^y^Aj^y^u^{x,t) = (3.13) 
= u^{x + yi'ti + 2/2 "^j, t) - u^{x - yi'ti + ?/2"e^j, t)- 
-u^{x + yi'ti - y2~tj, t) + u^{x + yi'ti + 2/2 "^j, t), 



where 



Denote also 



A.k,hu{x, t) = u{x + n't ki t) — u{x — h'tk, t). (3-14) 



r{,x\ t, yi, 2/2) = Ai^y^Aj^y^p{x', t). (3.15) 



Note that 



dx} dyf ' dx'j dyl 

Therefore in domains i?^ x {-00 < t < T} x {0 < yi < 1} x {0 < 2/2 < 1} the 
functions v'^{x,t,y) satisfy the equations 

. (3.17) 



d^v^ 3d^v^ 3d^v^ Id^v^ Id^v^ . 

+ T^r^ + T^r^ + T^r^ = 0, 



^. . 4 9x| 4 (9x2 '4 ^^2 4 



Note also, that 
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y2 



Exactly the same way 



and therefore 



where 



Similarly, we have 



Denote 



v^\<Ay^M{T), 



\v^\<^yrraM{T), 



Vmin = min{7/i,7/2} 



\r\< Ay^M{T), 



<Ay2N{T). 



(3.18) 



(3.19) 



(3.20) 



y^{yi,y2), ^ {{x,t,y) : \xm\ < l,m ^ 1, N - 1,0 < ±Xn < 1, (3.21) 



-1 < i < 0,0 < yjk < 1,A; = 1,2}. 



Denote also 



S± = dP^ \ {{t = 0} U {xn = 0}) , Eo = dP^ n {xn = 0} , (3.22) 



that is - are parabolic boundaries of the parallelepipeds P^ without their 
common part {x]\f = 0}, and the last will be denoted by Sq. 

In the parallelepipeds the functions and r satisfy the following problem 



L*v^^O, {x,t,y)eP^, 



(3.23) 



^^^|s± = g'^{x,t,y), 



(3.24) 



(3.25) 



n- el^x'T + 0^^ - — ^F{x,y,t), xn = 0, 



dx 



N 



dx 



N 



(3.26) 
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where 

F{x, t, y) = Ai,y, Aj- yj, g^{x, t, y) = Ai^y^Aj^y^u^, (3.27) 
and, in view of the assumptions (13. ip . 

\F{x,t,y)\<2{f)^]{2y,y. 
A similar inequahty with replacing j by i and ?/2 by yi gives similar to f l3.18p 



\F{x,t,y)\<CM{T)yl,,. (3.28) 



Note also that by flHTTSj) . 



\g^{x,t,y)\<CAf{T)y^,^. (3.29) 

To estimate and r, we're going to apply to the problem (I3.23P - (I3.26P the 
maximum principle in the following form. 

Lemma 3.3 Let functions H^{x,t,y) e C^'\P^)nC^''^(P^), S{x',t,y) e C^'H^o) 
satisfy the conditions 

L*H^>0, {x,t,y)EP^, (3.30) 
H^\j:± > 0, (3.31) 



+ A^S = 0, xn = 0, (3.32) 



Then 



St-eA^,S + h-^- r^<0, x^ = 0. (3.33) 

OXn UXn 



> 0, ix,t,y)eP ; S< 0, (x', t, y) G Sq. (3.34) 



We do not give a detailed proof of this lemma, since it uses standard argu- 
ments. We only note that the functions can not reach a negative minimum 
at {xat = 0}, as in this case, by f l3.32p . they would reached a negative minimum 
simultaneously and corresponding point would be, again by ( I3.32p . a point of 
a positive maximum of the function S. All this together in the standard way 
contradicts the boundary condition (I3.33p . 

We shall need the the auxiliary functions w^{x,t), defined on 

= {\xm\ < l,m = 1,A^- 1,0 < ±xn < 1,-1 < t < 0} 
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correspondingly, and such that 

dw 



.a , ^^^^ 3dV 3d'-w , , , , 

w^'^|{lxfc|=i}u{t=-i} > > 0, (3.36) 
y;±(0,0) = 0, w^|n±>0, (3.37) 

w^(x, t) G C^'^(n^ n {xtv = 0}). (3.38) 

Such functions can be constructed as follows. Consider for example, w^{x^t\ 
Let G+(a;,t) is a function from C°° in i?;!^ x (— c)0,c)o), such that = for 
\'^\ + N ^ 1/4 and for t < —2, |x| > 2 and > in the other points of 
X (—00,00). Let w;"'"(x, t) is the solution of the following initial boundary 
value problem in half-space 

L* = 0, xat > 0, t > -2, 
w^+|x.^=o = G+(a;,t)GC°°, 
tt^+|t=-2 = 0. 

Lemma 12.61 implies that the function exists in the appropriate class, and 

Because of the properties of G+(x,t) and by the strong maximum principle 
(see [M]), the function has all desirable properties, including fl3.36p . 
Now consider the following comparison functions defined on P . Denote 

<^(2/) = 2/1^/2 (2/1 + 2/2)"'^', (3.40) 
V'^(x^,?/) =7/17/2 [(1/1 ±a;7v)' + (z/2± a^Tv)']"^^'' ±a;jv>0, (3.41) 



1/min > 0, 



(3.42) 



The direct verification shows (cf. that the functions (p and ip"^ possess 

properties 

dip 



|a;iv=0 



0, 



dxN 

Thus, if we choose a sufficiently large constant > 0, the functions 
will have the properties 

L*h^ > i^lxNry;^^' > 0, {x,t,y) e P^, 

At the same time, the functions 6^{x,t,y) have the properties 

L*9^>0, {x,t,y)eP^, 



dxN 



mm; 
xjv=0 



^^lufc{|xfc|=l}U{t=-l} > ^Z/min- 

Consider now the following comparison functions 
H^{x, t, y) = [LiO^ix, t, y) + L^h^ixN, y)] {N{T) + M{T)), 

S{x', t,y) = - [L,9^{x', 0, t, y) + L^h^i^Q, y)] {N{T) + M(T)), 
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where Li h L2 are some positive constants. 

Choosing first Li and then L2 are sufficiently large, and using on one hand 
(ESHD, IK29\i . and on the other hand IK^ - (ESID, we see that the triple of the 
functions 

±v^{x,t,y), H~ ±v~{x,t,y), S±r{x',t,y) 
satisfies in to the conditions of the lemma ESI Hence, 

H+±v+>0, H-±v->0, S±r<0, 

that is 

\v^\<CH^, \r\<C\S\, {x,t,y)eP^. (3.57) 

Taking in fl3.57p x = 0, t = 0, in view of 6^{0,0,y) = 0, we obtain, for 
example, for similarly [28j 

I A,,, A,, ,,^+(0, 0)1 < C MT) + M{T)] y,y2 {y[ + y^) . 

Dividing both sides of this relation by yi and taking the limit with yi — )• 0, we 
obtain 



5^;(o + z/2^.,o)-^(o-y2V„o) 

and similarly 

^{0 + ya„0)-^{0-ya„0) 



<C[MiT) + MiT)]yl, (3.58) 



<C[Ar{T)+M{T)]yl (3.59) 



Since all of the above arguments are valid, as noted, for any (xq, to) ^ -R^~^'°, 
by the same token the estimate fl3.1ip and the lemma 13.21 are proved. 
■ 

We continue the proof of the theorem. It follows from (13. 5p . that 

pt - EA,,p = F{x\t) ^ -h^j— + + /, (3.60) 

Ox jv Ox TV 

Moreover, in view of the inequalities (ll.33p . (12. ip . (12. 4p 

< C (|V«+|;jL,.. + |Vu-|^|_,.) + CM{T) < 

<C{T + aYU{T) + CM{T), (3.61) 
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and 

F{x',0) = 0. 

Making in the problem (l3.6Up . (13. 6p the change of variables x' = £^^'^y, we 
obtain the problem 

Pt-Ayp = F{y,t), {y,t)eR^-''% p{y,-a) = 0, (3.62) 

and 



(3.63) 

It follows from the arguments of rji.IV, that 

w;:4L.<c ((%.,)-_. (%,,);-.,„). (3.65, 

Making in (I3.64p . (I3.65P the inverse change of variables, in view of 03.63P we 
obtain 



(/^*(^^^))!'<-- +^E(/^-»-'.>l';<-- ^ C\Fi^',t)\^iU.- (3.66) 
Thus, in view of the estimate (13.111) of the lemma 13. 2[ it is proved, that 



IPtl'^U. + {"^.'Py^^pL + ^\p\'^f-L < C{T + ayUiT) + CM{T), 
or, in view of (ll.20p . [26j and of the finiteness of p. 



IpI + e\p\ltX. < C{T + arU{T) + CM{T), (3.67) 

where the constant C does not depend on e > 0. 

Now, considering M^(x,t) as the solution of the Cauchy-Dirichlet problem 
(13. 3p . (13. 4p . (13. 6p . by the lemma [2l6] and the estimate (I3.67p . we conclude that 



uXnll^ + l^'l;,?:^ ^ + «)'^(^) + CMT). (3.68) 



I (2+7) , U,-|(2+7) 

It follows that in the condition (13. 5p 
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dx 



N 



C 



+ 



dx 



N 



C 



1+/3- 



< C{T+aYU{T)+CM{T). 



Thus, the function p{x',t) satisfies the Cauchy problem fl3.60p . (13. 6p with the 
right hand side F and the last has the property 



F{x', -a) = 0, „^ < C{T + afUiT) + CM{T). (3.69) 

Making again in (I3.60p . (13. 6 p the change of variables x' = s^^'^y, we arrive at 
the problem of the form (I3.62p with F, where the last is such that 

As above, completely similar to [23j, Ch.IV, for solutions of the problem (I3.62p 
we have the estimates 

{Pt)l^r-d <c{f) < C{T + aYU{T) + CM{T), (3.71) 

eV»j:;«T2<c(^)'';:::, (3.72) 

and we note that in obtaining the estimate (I3.7ip the condition F{x', —a) = is 
important. 

Proceeding as before and going back to the variables x', we find from (I3.72p 
and flHTnj) that 



N-l 

^ E {P^^^.)%P- ^ C (^)l!^f-t < C{T + arU{T) + CM{T). (3.73) 
Now combining the estimates (I3.73p . (I3.7ip . (I3.68P and (I3.67p . we find that 

l{{T)<C{T + a)''U{T) + CM{T). (3.74) 

Taking now in (I3.74p T = Tq, so that the value of To + a > is sufficiently 
small, we obtain estimate (13. 8p on the interval [— a, Tq]. Considering further the 
problem (13. 3p - (13. 7p on the interval [—a + (a + To)/2, —a + 3(a + Tq)/2] and 
removing the initial data with the known functions, that is moving along the axis 
of Ot up, exactly as in [23], Ch.IV, we obtain the assertion of the theorem 13. II on 
an arbitrary time interval [— a,T]. 

Thus, the theorem 13.11 is proved. □ 
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4 Reduction of the problem (ll.3l) -( frTf ) to the prob- 
lem in the fixed domain. 

Let p{uj,t) is the unknown function defined in Section [T] and parameterizing un- 
known (free) boundary Tp^T (p('^; 0) = 0), and let p{x,t) = Ep{u,t) is the exten- 
sion of this function to the whole domain Qt by the extension operator E from 

We pass in the problem (II. 3p - (II. 7p from the unknown functions u^{y,T) to 
the unknowns f ^(y, r) = \u^\^~^u^{y^ r). Then the relations (ll.3p - (ll.7p take the 
form: 

Lo{v^)v^ = ^ - a^lv^r'^yiy, r) = 0, {y, r) G fi^^, (4.1) 
v+iy,T) = v-iy,T) = 0, (y,r)Gr,,T, (4.2) 

N N 

a+^cos(l^,y,)v^^-a'^cos(l^,yi)v~ = kcos(l^,T), (y, r) G T^^t, (4.3) 
1=1 1=1 

v^{y, r) = \g^r-'gHy, r) ^ h^{y, r), {y, r) G r± (4.4) 

v^{y,0) = \u^r-'u^iy) ^ v^{y), y G TF. (4.5) 

We make in the problem (14. II) - (14. 5 p the change of variables (|/, t) = ep{x,t) 
which is defined in (11.311) . Denote for simplicity by the same symbols v^{x, t) the 
unknown functions after this change of variables, that is, 

v^{x,t) = v^{y,T) o ep{x,t). 

Then, in view of the properties of {y, r) = ep{x, t), in the variables (x, t) the prob- 
lem (I4.ip -( H3]) reduces to the following problem in the known fixed domains fl^ 
for the unknown functions v^,v~,p (besides x- variables we use the corresponding 
coordinates (w, A), which were introduced in (I1.28P ): 



Lp{v^)v^ = h^pPt - \v^rVlv^{x,t) = 0, (x,t) G (4.6) 



di 

v+{x, t) = v-{x, t) = 0, (x, t) G Ft, (4.7) 

^1 "^ii(a;,p)Pu;,P<.^,)(a^-^ - a"-^) = -^Pt(l + Pa), {x,t) ^Vt, (4.8) 

v^{x,t) = h^{x,t), (x,t)Gr^, (4.9) 

t;±(x,0) = t;^(x), X G n±, p(a;,0)=0, (4.10) 

p{x,t) = Ep{uj,t), (4.11) 
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where Vp = £pV and the matrix^^p is the conjugate and inverse to Jacobi matrix 
of the mapping (ll.3ip for t = const, mij{x,p) are some given smooth functions 
of their arguments, and 

Note that the last definition is legitimate, since the function p{x,t) is not identi- 
cally zero only if x G JV, where the coordinates (w. A) = {u^, Xx) of the point x 
are defined, and the coordinate A is independent of the choice of local coordinates 
u (we use the index {ux, Xx) to distinguish these coordinates for a point x from 
the corresponding coordinates {uy, Xy) for a point y). 

Below we explain the derivation of the relations f l4.6p -( BrT0|) . here we note the 
following. The relation ( 14. 8 p contains the expression 

N-l 

Sp = Sp{u,p,p^) = (1 + ^ mij{x,p)p^^p^^), 

which is explicitly expressed in the local coordinates u. But, in fact, the expres- 
sion Sp is strictly a function of the points of the surface Tt and its values at the 
points of Tj' does not depend on a choice of local coordinates u. Indeed, first, 
for any choice of local coordinates u the condition (14. 8 p is equivalent to (14. Sp . 
which is independent of a choice of local coordinates, and, secondly, all the other 
factors and the terms but 5*^ in the relation (14. 8p are invariant with respect to a 
choice of u and they are the function of the point of the surface only. Hence, 
the expression Sp, as a function of the point of the surface F^, is invariant on a 
choice of local coordinates u as well. And thus, the map p — ?■ Sp{u, p, p^)) defines 
a nonlinear operator, acting on functions defined on r^. This operator is invari- 
ant under choice of local coordinates uj, it acts in the space of functions on Tt 
and has a certain expression Sp{uj, p, p^J) for every particular choice of the local 
coordinates u. 

Further, the expression h'^Pt is the recalculated in the variables {x,t) 



derivative after the change of variables (ll.3ip : 

N-l 



dv"^ dv^ dt ^— ! dv^ dojxi dv"^ dXx 



dr dt dr ^ dujxi dr dXx dr 
1=1 



Here in fact 



dt duxi 

= ^ = 0. (4.13) 



and for the value of 4r^, due to the relation 



Xx = Xy- p{X,t) O Cp^, 
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and taking into account f l4.13p . we have 



dX d 

dp dt dp dXx yr-^ dp duxi dXr^ 



-Pt - Pa, 



_-|^ 'XI 



dt dr dXx dr ^ du^i dr dr 
So in the variables x and t 

dX^ 



Pt/(l + PAj. (4.14) 



dr 

Thus, it follows from fOB]) and fOij) . that 

dv^ .dv^ ., dv^ , , 

^ ° = ^ - [^/(l + p,)]p, = ^-K Pt. 

We explain further the transition from the condition (14.31) to the condition 
(14. 8 p under the change of variables f ll.3ip . as we shall need in the future the exact 
explicit form of this condition. Define in the neighborhood A/r of the surface Tt 
the function 



"^piy, r) = X^o ep^{y, t) = Xy - p{x, t) o ^(y, r) = X{y) - p{y, r), (4.15) 

where for simplicity we have retained for the function p(x, t) o ep^{y, r) the same 
notation p(y, r). By the definition ±$p(y, r) > for (y, r) E and $p(y, r) = 
for {y, t) G Tp^r- Hence in (Oj) 



cos(i^,y.) = ^^^, cos(i^,r)- ^''^ 



|V(y,r)$p| |V(,,.)$p| 

Therefore, the relation fl4.3p can be written as follows 

a+iVyV+, Vy^,) - a^iVyV-, Vy^p) = ■ (4.16) 

Under the change of variables f ll.3ip the right hand side of fl4.16p . due to the 
definition of $p, takes the form 

dX 

fcV = ^^ = -W(l + PAj, (4.17) 

owing to f l4.14p . 

On the other hand, under the change of variables (ll.3ip 

(Vyv^, Vy^p) o ep{x, t) = {Vpv^, VpA,). (4.18) 

Denote by A(x) the transition matrix from the gradient with respect to the vari- 
ables X to the gradient with respect to variables {u^, X^), that is 
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V, = A(a;)V(A.,^,) (V, = A(y)V(A,„ 



UJy)), 



where 



A(a;) 



/ dX duji 
dxi dxi 



dxi 



dX dug du)M-i 
\ dxisi 9a: jv "' <9a;jv / 



(4.19) 



(4.20) 



and similarly for the variables y. Then in the variables (x, t) 



Note that V (\^^^^)Xx = {1,0,...,0}, and also t>='= = on F, hence dv^/dui = 0, 
and therefore 



± r n n1 ri n n1 V7 



{1,0, ...,0} 



5A, 



{Xx,uJx)Xx- 



Thus we obtain 



{Vyv^, Vy^p) o ep(x, t) = {Vpv^, VpA,) = ^(VpA„ VpA,). (4.21) 
On the other hand, due to the definition of $p(?/, r), 

(VpA,, VpA,) = (Vj,(A, o e;^), Vj,(A, o e^^)) o Cp = (Vj,$p, Vj,$p) o Cp. (4.22) 
Using introduced in (I4.19P matrix A{y), we have 



(4.23) 



(V,$p, V,<fp) = (A(|/)V(A„^^)<l'p,A(y)V(A„a.,)$p) = 
= (V(A„c.,)$p, A(|/)*A(y)V(A„.;,)$p). 

First, by the definition of $p, 

9a; = 9a;^^^-^^^'")) = '-^^- 

In addition, since the coordinate Xy is counted by the normal to F, and Uyi are 
coordinates on the surface F, then 



(4.24) 



(V,A(y), Vy\{y)) = 1, (Vj,A(y), V,u;,(y)) = 0, z = 1, iV - 1. 
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Therefore the matrix A*{y)A{y) has the form 
A*(y)A(y) 



/ 1 

mil 





mi2 



\ 

"^l(Af-l) 



(4.25) 



\ m(jv_i)i m(Ar_i)2 ... m(Ar_i)(7V-l) / 

where 

niij = rriji = (VyUiiy), VyUj{y))- (4.26) 

are some smooth functions. 
Thus, 

(V(A„.,)$p, A*(y)A(7/)V(A„.;,)$p) = 

N-l 

= (1 - PaJ' + Yl ^^3iy)P^y.P^yr (4-27) 

Make now in fl4.27p the change of variables f ll.3ip . and recalculate the deriva- 
tives of p with respect to {Xy^Uy) in terms of the derivatives with respect to 
(Aa;, ojx)- We have 



dt_ 



N-l 



^'^y i=i 

It follows from the definition of the mapping Cp that 



dt 



0, 



9X 



0. 



(4.28) 



(4.29) 



At the same time by (g^HD, 



that is 



1 - Pa = 1 - Pa, 



OK 



d\y 1 + Pa, 



Therefore by fOH]) . fICTD and (ICTj) 

Pa, ° ep 

Further, 



Pa, 



1 + Pa, 



(4.30) 



(4.31) 



dt 



Pujy, oep = Pt 



N-l 



dio. 



(4.32) 



yi 
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and 



dt 



yi 



0, T^ = 5.„ 2,J = l,...,iV-l. 



At the same time 

d{\x o ep) 



duj. 



yi 



d 



(Ay-p(i/,r)) 



That is by virtue of and 
hence by (Olj) . 



Pec 



Thus, it follows from f lOT]) . (H^ZD, f lOT]) and that in flOT]) 

(VpAa;, VpAi 



.1 + Pa, 



1 + ^ Tnij{x,p)p 



(4.33) 



(4.34) 
(4.35) 

(4.36) 



Finally, the relation follows from the relations (HTTB]) . (HTT]) . and 



5 The linearization of the problem. 

Our goal in this section is the extraction of the principal linear part of the problem 
fl4.6p - fl4.lip in terms of the deviation of the unknown functions (f"^,f~,p) from 
functions constructed from the initial data and satisfying fl4.6p - fl4.lip for t = 0, 
as it was done in [3j, [6J. 

Note that from the equations (14. 6p . (14. 8 p and from the initial data f l4.10p we 
can calculate the derivatives with respect to time dv^ jdt and dp/dt at t = 0: 



dp 1 
-(c.,0) = p,(c.)^-(a 



dX 



^(x,0) = vHx) ^ ^p, + a± |.±(x)r VV(^), 



(5.1) 
(5.2) 



and, in view of the assumptions fll.lOp . f ll.26p . 

p^{co) e Ci+^'-"(r),/3' ^ y(l - «/2) > /3, vtix) G C]'(JW). (5.3) 



Completely analogous to Ch.IV, on the base of results Pl] on the solv- 
ability of the Cauchy-Dirichlet problem for degenerate equations we construct 
such functions w^{x,t) G Cs^"' ^"^{^l^) that 
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and 



w^{x, 0) = v^{x), (x, 0) = t)|rT = 0; t)|p± = /i^- (5.5) 

In addition, just as described in [23j, Chapter 4, there is a such function a{u), t) G 
p3+„._<,.,+l±tajp^j that 

kir-"'"*^'<c(i%irji?''+i%irjp'') (5.6) 

and 

a{u, 0) = p{u, 0) = 0, ^(a;, 0) = pi{u). (5.7) 

Moreover, by the method described in [23j, Ch.IV, the function a{u),t) can be 
extended with the class and with the inequahty (15.61) to a function defined in Qt 
which is non-zero only in the neighborhood Af x [0,T] of the surface Ft- 

The linearization of the relations f l4.6p - fl4.lip consists in the following (we 
describe the general scheme of the arguments - the exact formulations will be 
given below). We denote the space 



p2+/3-"(r^) = {p:pe c'+^-^-'-^iTT), pt e ^^+^-'^(1^)} (5.8) 

with the norm 

|p|p2+,-.(r^) = |p|^,+,_,2+^^^^^ + (5.9) 
Denote also 



V'O = (u7+,U7~,(t), 

and represent the the relations fl4.6p - fl4.lip as 

F{i)) = Q (5.11) 

with some non-linear operator of ip. Keeping essentially in mind the application 
of Newton's method, we represent the relation (15.111) as 
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= fo + G{^-tPo), (5.12) 

where F'{ipQ) is the Frechet derivative of F{ip) at the point ipQ. In this 
the new unknown we consider the difference 

if = ip — ipQ = iy^ — — , p — a)^ (5.13) 

which belongs to the spaces with zero, that is, 

^en,^ <r''''"/'(^?) X cZ'^'^''\n=^) X p',^<'-\Tt). (5.14) 

By the construction of the element ipQ = {w~^, w~, a), it has an increased smooth- 
ness (7' > 7) and satisfies the relation -F(V'o) = for t = 0. Therefore, using the 
inequalities fl2.ip - fl2.4p . we can estimate 

ll/oll = ||-F(^o)|| <C7T^ (5.15) 

Below we show that the operator F'{iIjq) has the bounded inverse in the ap- 
propriate spaces, so that the equation f l5.12p can be rewritten as 

ip=[F'{^o)]-'fo + [F'{M-'Giv)^ 

= ho + H{ip) = K{ip), (5.16) 

where by (I5.15P 

\\ho\\ < CT^", (5.17) 

and the operator H{(f) is the "quadratic" with respect to (p by the smoothness 
of F{iIj) in its argument and by the definition of G{ip — ipo) = G{(f) in (I5.12p : 

\\H{y,)\\ < C , \\H{ip2) - H{ip{)\\ < C{\M\ + 11^211) \\V2 - ViW ■ (5.18) 

For sufficiently small T > it follows from f l5.17p and fl5.18p that the operator 
K{ip) maps some small ball Br C with a small r into itself and K{(p) is 
a contractive there. The only fixed point of this operator gives, obviously, the 
solution of the original problem. 

Thus, our goal now is to write the problem fl4.6p -( HTT]) as fl5.12p . 

Denote 

e^ = v^-w^, 6 = p-cT. (5.19) 
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Lemma 5.1 The problem ( I4.6p - fl4.1ip can be represented as a problem for the 
unknown functions 9^ and 5 as follows 

= (x, t; 9, 6) + F^{x, t- 9, 6) = F^{x, t; y?) + F^{x, t; (x, t) G Q^, (5.20) 



^+ = ^+ = 0, {x,t)eTT, (5.21) 



= F3{x,t-^) + F4{x,t;^), {x,t)eTT, (5.22) 

^^ = 0, {x,t)eT$, (5.23) 

9^{x,0) = 0, S{uj,0) = 0, (5.24) 

5{x,t) = E5{uj,t), (5.25) 

where for arbitrary (p = {9~^,9~,6) & Br C Ho, r < 7o/2 in the righthand sides 

Fi of the relations fl5.20p - fl5.24p all the functions Fi vanish at t = and the 
following estimates are valid 

\Ft{x,t;^)\^jl^<CT^, (5.26) 

\Ftix, t- ^2) - F±(x, t- ¥.01^"^ < CT^ \W2 - , (5.27) 

|F2±(x,t;^)|(|^<Cl|v.l|^, (5.28) 

\Ft{x,t- ^2) - F^{x,t- ^i)|(|^ < C{\MW + WviWu) \\V2 - ViWn ' (5-29) 

|F3(x,t;^)i;.';^-"''^^ < CT^ (5.30) 

|F3(a;, t- ^2) - F3±(x, t; ^OI^'J""""''^^ < CT'^ ||^2 - , (5.31) 

44 



+ m(1+'5-".^^^) ^ ^11 ||2 



(5.32) 



Proof. 

Meaning of the inequalities (I5.26p -( 15.33p is that, according to (I5.12p . the ex- 
pressions and F3 contain smoother terms and to evaluate them, we use in- 
equalities fi2.ip - fi2.4p . and the expressions and F4 are "quadratic" with respect 
to (p. 

In the case of a uniformly parabolic equation in f l4.ip (f lTT^ ). this lemma is 
proved in details in [B], Section 2.3. Therefore, we mention only the differences 
that arise in the case of degenerate equations. 

First, in contrast to the [6], we can not expect that the extended function 
p{x, t) = Ep{uj, t) satisfies the condition dp{x, t)/dX = on Ft, and so we explain 
the obtaining of the relation fl5.22p from the relation f l4.8p . The relation fl5.22p is 
obtained from f l4.8p explicitly after substitution in fl4.8p the expressions = 6^, 
p = 6 + a and the transfer of the junior and quadratic terms in the righthand 
part. It is easy to verify that (15.221) coincides with (14.81) for 



F3{x,t,(p) 



'N-l 



dv^ _ dv 



kaf — (a 



~~dX 



_ dw 
~~dX 



k(Jt — {0. 



~~dX 



_ dw 
~~dX 



[kpt(Tx]} 



(5.34) 



F4{x,t,ip) = -k6t6x, (5.35) 

where in fl5.34p p = 6 + a, v"^ = 0^. 

Using the fact that 5{uj, 0) = o"(w, 0) = p{uj, 0) = 0, estimating each term in 
square brackets in fl5.34p separately, and using the inequalities fl2.ip - fl2.4p it is 
easy to obtain for F^{x^t,ip) the estimates fl5.30p . f l5.3ip . For example, since by 
the construction 



/ ,dw^{x,Q) dw (a;,0), ^ 
kat{x,Q)-{a+ wr^-oT ^^-^) = 0, a; G F, 



dX 



dX 



then 
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kat 



~~dX 



_ dw 



(l+/3-a,l±^) 



< 



kat — (ct 



~~dX 



_ dw 
~~dX 



l+/3-a \ 

/ 



< 



Since this term is linear with respect to ip, this yields fl5.30p . fl5.3ip for this term. 
The remaining terms in the definition of F^{x,t,(p) are treated similarly. 

As for F4^{x,t,ip), the estimates fl5.32p . fl5.33p for this expression are obvious 
because it is quadratic. 

Another difference from [6j is the presence of a degenerate factor in the third 
(elliptical) terms in the left-hand side of (14. 6p . Represent this term as (we consider 
only the equation for the sign " + ") 

{v+rVlv+ = «)°V>+ + [(t;+)" - «)lV2t;+ = 



— V^O / " p 

and show that A^l^p) satisfies the inequalities 



(7) 



(5.36) 



(5.37) 



Let us assume that the function A = A(x) is extended with the preservation 
of the class from the neighborhood A/" of the surface F on all Q to a function 
satisfying the conditions 

ud+{x) < X{x) < u-^d^{x), 
retaining for her the same notation. Write further ^"'"((y9)in the form 



Since vq, f = for x G T, then for x G A/" 



A"(x)V;w^ 



(5.38) 



s, u)ds, 



(5.39) 



where we assume the neighborhood TV so small that 
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dX 



^ > z/ > 0, X eX, 



(5.40) 



In addition, we assume that T is so small that 



dw 



— ^x,t)>u>0, {x,t) e Af X [0,T]. 



Assuming now that the radius r = r(z/) of the ball Br 3 is sufficiently small, 
we can assume that for if G Br 



^{x,t) = ^{x,t) + ^{x,t)>u>0, {x,t)eNx[0,T]. (5.41) 

In addition, outside the neighborhood Afr holds 

v^{x)>u>0, xeIl+\Af. (5.42) 
Therefore, assuming as above T and r sufficiently small, we can assume that 



v+{x, t) = w+{x, t) + ^+(a;, t) > z/ > 0, x eVL+\ N. 
Thus, we have the representation 



(5.43) 



f ^(x, t) 



and 



and also 
Similarly 



so 



l$+(x,t,^)l^;-^<c( 



dv' 



(7) 



^^{x,t,(p) > z/ > 



X{x) 



$+(x) 



"0 

A ' 



xen+\Af, 



<^+ix] 



< C, <^^{x) >u>0. 



In addition, it follows from the representations (I5.44p and (I5.47P that 

|$+(x,t,v.)-$+(a;)|^^lp< 



(5.44) 

(5.45) 
(5.46) 

(5.47) 
(5.48) 
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(7) 



<Ciy\\)T>^, (5.49) 



and similarly 



|<l>+(a;,t,(^i) - <|.+ (x,t,(^2)|^|T < C^'' 11^2 - ViW . (5.50) 



By the properties fl5.46p and fl5.48p . the mapping 
is smooth, and by (I5.49p . (I5.50p 

|5i(y.2) - 5i(^i)|1|t < 11^2 -^iWn- (5.52) 

Now from the definition of the expression A'^{ip) in fl5.38p .from the relations 
(I5.52P and from the smooth dependence of Vp = £pV on p, taking into account 
that the factor A"(x) is appropriate for the weighted estimates of the second 
derivatives of and in the space Cq~^"''^~^^'''^{^1'^), it is easy to see that 

\A+{y^2) - ^'"(</'i)|1|t < 11^2 - ViWn ■ (5.53) 

Thus, in view of fl5.36p and f l5.53p the linearization of the equation (14. 6 P is 
reduced to the linearization of a linear on equation that was done in details 
in [0]. 

Note also that the insignificant difference between (I5.20p from [6] is still in 
that we, in fact, leave in the left-hand side of (I5.20p only the leading terms, 
moving all the other to the expression F^(x, t, ip). 

This completes the proof. ■ 



6 The linear problem corresponding to the prob- 
lem (Km-iKm. 



In this section we consider the linear problem obtained from the problem (I5.20p - 
(I5.24P for a given right-hand sides from corresponding classes. In this case, \v^{x) \ 
is replaced by d^(x)B^{x,t) ~ X{x)d\V^{x). And, as in the previous section, we 
assume that X{x) is extended to all to a smooth function of the class H^^"' , 
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u < X{x),d^{x) <u-\ xen\Af. (6.1) 

Thus, we consider in the domains fi^ the following problem of finding the 
functions v^{x,t), defined in the domains fi^, and the function 5(0;, t), defined 
on Ft, on the conditions 



^""^ - A(a;)"5±(x,t)V2i;± - A^{x,t) - X{x)''B^{x,t)V^6 



dt ' ' ' ' ' ' ' ' \dt 

= fHx,t), {x,t)en$, (6.2) 

v+{x, t) = v'{x, t) = 0, (x, t) G Tt, (6.3) 



Ov~^ dv 
k6t — eArS + ( a'^—, 



dX 



dX 



6x{a^A+{x,t) -a'A'{x,t)) 



/2(x,t), {x,t)eTT, 



v'^{x,t)=0, (x,t)er 



v^{x, 0) = 0, 6{x, 0) = 0, X e 



(6.4) 
(6.5) 
(6.6) 



6{x,t) = E6{u,t), (6.7) 

where the extension operator E was defined in the section O Ap is the Laplace- 
Beltrami operator on the surface F ( compare flj). We assume that 

fHx,t) G Q;/'(fi^), f2{x,t) G Co^"^" '^(Ft), (6.8) 
e, a^, k are given positive constants, 

u <k, a^, B^{x, t), A^{x, t) < z/-\ (6.9) 

A^{x,t),B^{x,t) G C^'^'^iS)^). (6.10) 

For the problem fl6.2p - fl6.7p by the standard method of the freezing of coeffi- 
cients and multiplication by smooth cutting functions we can obtain the Schauder 
a priori estimates of the solution completely similar to |23) (or [6] in the case of 
the Stefan problem). At that the model problem, obtained by the freezing of 
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the coefficients in points of the boundary F at t = 0,with the subsequent local 
rectification of the boundary, was studied in the section [31 At considering such a 
model problem the functions B^{x,t) and A'^^xt) are replaced by the constants 
B"^ = B^{xo,0) and A"^ = A^{xo,0), xq G T. After this, the change of the 
unknown function 

u^{x,t) = v^{x,t) - A^5 (6.11) 

reduces the problem (16. 21) - (16. 7p with the frozen coefficients and with the flat 
boundary exactly to the problem (I3.3l) - (l3.7p . 

From these model problems associated with the boundary F we get the es- 
timate of the function 5(x, t)|r and border estimates of the functions v^{x,t). 
After that, the rest of the model problems associated with a strictly interior 
points of are standard because of the condition (16. 7p . and due to the absence 
of degeneracy of the equations at these points - see |23j. 

Thus, the following is true. 

Lemma 6.1 Suppose that the conditions (I6.8p - (l6.10p . Then for the solution of 

the problem (E2])-(EZD from the class E Cs~^^'^ (U^), 6 G C3+^-"'1+^^^(Ft) 
the following a priori estimate is valid 



N—l 



< Ct (\ft\%. + |/r|S- + Ihtr:'"""^') - CrMin (6.12) 

where the constant Ct in (I6.12p does not depend on e E (0, 1). 

We now show the solvability of the problem (I6.2p - (l6.7p . 

Theorem 6.2 Suppose that the conditions (I6.8p - (l6.10p are satisfied. Then for 
e G (0,1) the problem (16. 2p - (16. 7p is solvable in the space E Cs ' ^ (^t); 

5 G (73+/3-a,n 2-Q {Vt), o,nd the estimate of the solution (I6.12p is valid. 

/iTi-mv /iTi-^ -I- 2+7, — ± 

When e = the problem (I6.2p - (l6.7p is solvable in the space v E Cs ' ^ {^t)> 

6 G C2+''-"'^^^(Ft) a G C^+''-"'^^(Fr), and the estimate (KI^ without 
the term with e is valid. 

Proof. 

Define the linear operator M : 5 — )■ f ^ — )• MS which maps a function 5 G 
^2+/3-a, (Tt) fist to the functions f^, as the solution of the problem (16. 2p . 
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(16. 3p . fl6.5p . (16. 6 p with the given function S in fl6.2p . and then the functions the 
operator M maps to the function MS, which is determined from the condition 
(16. 4p with the given and 6x, that is the function M6 is the solution of the 
problem 

a^-g^ - ^'-g^) + ^A(a+A+ - a-A-), (6.13) 
M5{uj,0) = 0. 

By the theorem 12. TI this operator is well defined, and with £ > 0, by the known 
properties of the problem (I6.13P , 

|^^|,3+,-„,,.^, ^ + M(T)) , (6.14) 

- M^ir-"'""*^' < C.^T Ife - ^.ir"'"*^' • (6.15) 
Consequently, by (12. 4p . 

<C7,,rT^|52-5i|r, ' (6.16) 

Thus, for a sufficiently small T = the operator M is a contraction on 

Cg ' (Ft) and therefore has a unique fixed point, which by (16.141) . belongs 

also to the space Cq ' (F^) and together with the corresponding 

gives the solution of the problem. The estimate of the solution is given by the 
lemma [6Tl Moving now step by step up the axis Ot as in |23], we obtain the 
theorem with e > for any T > 0. 

Further, by the estimate (I6.12p . considering the sequence of the solutions vf, 

(^g, £ — )■ 0, we see that this sequence is compact in the spaces Cq^, ' ^ (^t) 

2 I ^ a, ^+^~" _ — _ 

Cg ' (Ft) correspondingly for any 7 < 7, /3 = 7(1 — a/2). The passing to 

2+-, 2+3 —± 2+;a-a,^+^~" 

the limit of this sequence in the spaces Cg^ ' ^ (^t) C*g '2a (p^^^ ^ 

Cg ' (Ft) gives the solution of the problem (16. 2p - (16. 7p for e = 0. Besides, 
as it follows from the uniform in e estimate (I6.12p . the limit function and 5 

belong to the spaces Cg't ^ (1]^) and Cg ' (Ft), 5t G Cg ' '"'^ (Ft) 
cor r esp ondingly. 

Thus, the theorem 16.21 is proved. 
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7 Completion of the proof of the theorem 



We complete the proof of the theorem 11.11 according to the scheme described in 
the section [51 

Define on the space "Ho (^o is defined in (I5.14p ) a nonhnear operator J^{if), 
if = {6~^,d^,6) G Hq in (I5.20p - (15. 25 p . which maps a given (f in the righthand 
sides of the nonhnear relations ( I5.20p . (I5.22p to the solution of the linear problem 
defined by the left-hand sides of these relations. It follows from the theorem 16.21 
and the lemma 15. 1[ that the operator J^{(p) has the following properties on the 
ball Br = '■ llv^llw ^ C Ho with the sufficiently small radius: 



It is easy to see that It follows from the relations (17. ip and (17. 2p that for 
sufficiently small T and r the operator J^(v?) maps the closed ball into itself 
and is a contraction there. The only fixed point of this operator gives the solution 
of the original nonlinear problem with free boundary. Thus, the theorem 11.11 is 
proved. □ 



\\Hv)\\n<C{T^ + T)\MH + CT^, 



(7.1) 



||^(<^i) - J-(<^2)lk < C{T^ + r)\\vi-V2\\H. 



(7.2) 
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